Accepted January 16, 2012 

Preprint typeset using I^TfTpC style emulateapj v. 5/2/11 



EVOLUTION OF MAGNETIC FIELD AND ENERGY IN A MAJOR ERUPTIVE ACTIVE REGION BASED ON 

SDO/RMI OBSERVATION 

Xudong Sun 1 ' 2 , J. Todd Hoeksema 1 , Yang Liu 1 , Thomas Wiegelmann 3 , Keiji Hayashi 1 , Qingrong Chen 2 , Julia 

Thalmann 3 

Accepted January 16, 2012 

ABSTRACT 

We report the evolution of magnetic field and its energy in NO A A active region 11158 over 5 days 
based on a vector magnetogram series from the Helioseismic and Magnetic Imager (HMI) on board 
the Solar Dynamic Observatory (SDO). Fast flux emergence and strong shearing motion led to a 
quadrupolar sunspot complex that produced several major eruptions, including the first X-class flare 
of Solar Cycle 24. Extrapolated non-linear force-free coronal fields show substantial electric current 
and free energy increase during early flux emergence near a low-lying sigmoidal filament with sheared 
kilogauss field in the filament channel. The computed magnetic free energy reaches a maximum of 
~2.6 x 10 32 erg, about 50% of which is stored below 6 Mm. It decreases by ~0.3 x 10 32 erg within 
1 hour of the X-class flare, which is likely an underestimation of the actual energy loss. During the 
flare, the photospheric field changed rapidly: horizontal field was enhanced by 28% in the core region, 
becoming more inclined and more parallel to the polarity inversion line. Such change is consistent with 
the conjectured coronal field "implosion", and is supported by the coronal loop retraction observed 
by the Atmospheric Imaging Assembly (AIA). The extrapolated field becomes more "compact" after 
the flare, with shorter loops in the core region, probably because of reconnection. The coronal field 
becomes slightly more sheared in the lowest layer, relaxes faster with height, and is overall less 
energetic. 

Subject headings: Sun: activity — Sun: photosphere — Sun: corona — Sun: surface magnetism 



1. INTRODUCTION 

E xtreme solar a ctivity is powered by magnetic energy 
(e.g. lForbesll200Q[ ). Flux emergence and shearing motion 
introduce strong electric currents and inject energy to 
the active region (AR) corona. Coronal fields are sub- 
sequently reconfigured, accumulating large amounts of 
magnetic free energy. When the overly energetic field 
gets destabilized, part of its excess energy is released 
rapidly, with power enough to drive explosive phenomena 
such as flares and coronal mass ejections (CMEs). Dur- 
ing this process, fast, irreversible changes of the photo- 
spheric field have been observed as the possible imprint of 
the coronal activity, including transverse field and mag- 
netic shear increase (e.g. Wang et al. 1994; Wang 2006) 
and longitudin al field decrease ([Sudol fc Harvevl 120051 ). 
Hudson! (|200Q[ ) suggested that the energy loss during the 
explosion causes "implosion" of the co ronal field, which 
is supported by r ecent observations (|Ji fc Wangl 120071 : 
iLiu fc Wandl201(l . 

This "storage-release" picture provides a scenario suc- 
cessful in exp laining many observed phenomena (e.g. 
iSchrijverl 120091 ). However, detailed knowledge of the 
process is still lacking, especially in a quantitat ive spa- 
tially and temporally resolved manner (e.g. iHudsonl 
|2011[ ). In this sense, uninterrupted, frequent photo- 
spheric vector field observation may be crucial, as it di- 
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rectly monitors the AR's evolution and provides infor- 
matio n on its non-potential nature (e.g. iSchrijver et al.l 
2005). In addition, field extrapolation models based on 
the vector boundary may provide valuable diagnostics 
of th e changing coronal field (e.g. [B egnier & Ca nfieldl 
2006; iThalmann fc Wiegelmannl 120081 : Uing et al.l l2QQ9h 7 
The re cently launched He liospheric and Magnetic Imager 
(HMI; iSchou et al.l |2011[ ) on board the Solar Dynamic 
Observatory (SDO) therefore presents a unique oppor- 
tunity to better understand the AR energetics with its 
full disk, high resolution and high cadence vector mag- 
netograms. 

Here we report the evolution of magnetic field and en- 
ergy in NOAA AR 11158 using a series of HMI vector 
magnetograms and a non-linear force free field (NLFFF) 
extrapolation. This 5-day uninterrupted, 12-minute ca- 
dence data set allows us to study in detail both the long- 
term, gradual evolution, as well as the rapid changes dur- 
ing an X-class flare and CME eruption. We briefly de- 
scribe the data set and the extrapolation method in Sec- 
tion O The long-term evolution and the rapid changes 
are discussed in Sections [3] and 2J respectively. In addi- 
tion to the HMI data, we employ coronal and chromo- 
sphere images, such as t hose taken by the Atmospheric 
Imaging Assembly (AIA; iLemen et al.ll2011[ ) on SDO for 
context and validation of the results. 

2. VECTOR MAGNETOGRAMS AND FIELD 
EXTRAPOLATION 

The HMI instrument on SDO observes the full Sun at 
6 different wavelengths and 6 polarization states in the 
Fe I 6173 A absorption line. Filtergrams with 0.5" pixels 
are collected and converted to observable quantities on a 
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Fig. 1. — Observations and modeling results for 2011 February 14 20:35 UT, about 5 hours before the X-class flare, (a) Remapped HMI 
vector magnetogram for the center region of AR 11158 as viewed from overhead. Vertical field (B z ) is plotted in the background; blue 
(red) arrows indicate horizontal field (Bh) with positive (negative) vertical component. Contours are plotted at ±600 G. (b) Vertical 
current density (J z ) derived from 5-pixel Gaussian-smoothed vector magnetogram. Contours are for B z and are identical to (a). Part of 
the polarity inversion line (PIL) is plotted as thick cyan curve, (c) Image from AIA 171 A band showing the corona magnetic structures. 
The dotted box in the center indicates the FOV of Figure [2Jl. (d) Selected field lines from the NLFFF extrapolation plotted over a cutout 
from the vertical field map. The lines are color-coded by the vertical current density at their footpoints (see the color bar); red field lines 
correspond to strong current density. The white dotted box indicate the FOV of (a) and (b). The FOVs of (c) and (d) are identical, about 
218 x 218 Mm 2 , or 302" x 302". Features of interest are marked in each panel; see text for details. (An animation of the vector field data 
set for the entire 5-day period is available in the online journal.) 



rapid time cadence. For the vector magnetic data stream, 
each set of filtergrams takes 135 seconds to complete. 

To obtain vector magnetograms, Stokes parameters are 
first derived from filtergrams observed over a 12-minute 
interval and then inverted through a Milne-Eddington 
based algorithm, the Very Fast Inver sion of the Stokes 
Vector (VFISV; iBorrero et al] l20TTh . The 180° az- 
imuthal ambiguity in the transverse field is resolved by 
an improved version of the "minimum energy" algo rithm 
(|MetcaHll99llMetcalf et al.ll2QQ6HLeka et al.l [2Q09). Re- 
gions of interest (ROIs) with strong mag netic field are 
autom atically identified near real time (jTurmon et al.l 
l2010f ). A detailed description on the production and rel- 
evant characteristics, spec ifically for the d a ta set for AR 
11158, can be found in iHoeksema et al.l (j2012[ ). Out- 
standing questions regarding data reduction are also de- 
scribed there. 

The magnetic field above the lower chromosphere 
large ly satisfies the force-free criteria (|Metcalf et al.l 
Il995f h where currents align with the magnetic field. For 



the most general case, 

VxB = aB, (1) 

B • Va = 0, (2) 

with the torsion parameter, a, varying in space but con- 
stant along each fie ld line. We use an optimization-based 
NLFFF algorithm (Wiegelmann 2004) to extrapolate the 
coronal field from the magnetograms in a Cartesian do- 
main. A preprocessing procedure removes most of the 
net force and torque from the data so the boundary 
can be more consistent w ith the force-free assumption 
(j Wiegelmann et aH 12006) . For reference, we also con- 
struct a potential field (PF) from the same observation 
using the vertical compo nent of the fi eld and a Green 
function algorithm (e.g. Sakurai 1989]). The magnetic 
free energy (Ef) can be inferred by subtracting the PF 
energy (Ep) from the NLFFF energy (En), where the 
energy is computed from the field strength within a cer- 
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tain volume V: 

E f 3_ dV _ f 3i dVm (3) 
j Jv ^ Jv ^ 

Here the subscripts N and P denote NLFFF and PF 
respectively. We note that the VFISV inversion scheme 
for HMI has the magnetic filling factor held at unity, so 
the obtained "magnetic field" is essentially averaged flux 
density. In the context of this work, we do not distinguish 
these two and use the unit Gauss (G) throughout. 

We use 600 vector magnetograms of NOAA AR 
11158 with a cadence of 12 min, from 2011 Febru- 
ary 12 to 16 (see online animation of Fig. [1] for the 
data set). The images are de-rotated to the disk cen- 
ter and remapped using Lambert equal area projection 
(jCalabretta fc Greisenll200l lThompsonll2006D . The field 
vectors are transformed to Heliographic coordinate s with 
projection effect removed (Gary & Hagyard 1990). For 
direct analysis of the photospheric field, we use full res- 
olution data at about 360 kmpix -1 (0.5"). For extrapo- 
lation, we bin the data to 720 kmpix -1 (about 1") and 
adopt a computation domain of 216 x 216 x 184 Mm 3 
(300 x 300 x 256). 

Because the preprocessed bottom boundary emulates 
the magnetic field in the chromosphere, we assign a uni- 
form altitude of 720 km (1 pixel) to thi s layer where 
the f ield has become largely force-free (|Metcalf et al.l 
1995) and use it through out the study. In general, we 
limit the field-of-view (FOV) for analysis to the center 
184 x 144 x 115 Mm 3 , covering most of the strong field 
region. Uncertainties are reported as mean ± standard 
deviation unless specified otherwise. 

The extrapolation is performed at 1-h cadence. For 
the 9 h around the X-class flare the full 12-min cadence 
is utilized. For convenience, we choose February 15 00:00 
UT as time (T = h) and use this convention when 
needed. The AR passed central Meridian on early Febru- 
ary 14 (-22.2 h). The X2.2 flare (W21S21) started at 
1.7 h (February 15 01:44 UT) and peaked at 1.9 h (01:56 
UT) in the GOES soft X-ray (SX R) flux and was accom - 
panied by a front-side halo CME (Sch rijver et al.l l2011). 

3. LONG-TERM EVOLUTION 

3.1. An Illustrative Snapshot 

We use a snapshot taken at 20:35 UT on Febru- 
ary 14, about 5 h before the X-class flare, to illus- 
trate the magnetic field structures in AR 11158 when 
it is well developed. It mainly consists of two bipoles 
(P0/N0 and Pl/Nl in Fig. QJi, P denoting positive and 
N negative) that form a complex quadrupolar structure 
(jSchriiver et al.|[20TT[ ). The total unsigned magnetic flux 
(|<£|) is about 2.7 x 10 22 Mx and the flux is balanced to 
within 1%. Vertical field (B z ) in the center of umbra can 
be as strong as 2600 G. Strong shearing motion exists 
between PI /NO as well as amongst a few newly emerged 
pores. The horizontal field (Bh) in the AR center is 
largely parallel to the major polarity inversion line (PIL; 
Fig.QJi), near which strong vertical electric currents (J z ) 
are present (Fig. [TJd). On the south side of the PIL, an 
elongated strip of positive flux following PI forms a so- 

4 http: //sun. Stanford. edu/~xudong/Article/field.mp4 



TABLE 1 

Estimated field parameters in the filament 

CHANNEL. 



Parameter^ 


Unit 


Cross section b 


Near axis c 


A 


Mm 2 


59.9 


4.2 




10 20 Mx 


4.87 


0.49 


(B) 


10 3 G 


0.97 ±0.27 


1.20 ±0.05 


(M> 


degree 


26 ±16 


11 ±8 


(Jh) 


mA m -2 


19 ±9 


43 ±5 


(a) 


Mm" 1 


0.27 ±0.15 


0.46 ±0.05 


(L) 


Mm 


29 ±5 


24 ±3 




Turns 


0.64 ±0.17 


0.92 ±0.17 



a All parameters are estimated from uniformly sampled 
points on a vertical cross section through the NLFFF ex- 
trapolation domain, marked as "CS1" in Figure [2] Angle 
brackets refers to the mean value. The notations are as 
follows: cross section area A, axial flux <E> a , field strength 
B, inclination angle 7, horizontal current density J^, tor- 
sional parameter a, loop length L, twist angle (j>. The 
inclination angle is measured with regard to the photo- 
sphere, ranging from —90° to 90° with 0° for horizontal 
field and sign consist ent with polarity of B z . a is esti- 
mated using equation lBll <f> is divided by 2tt to show the 
estimated number of turns. 

b Values are calculated within the Jh > 10 mAm -2 
contour in Figure [2^ and reported as mean ± standard 
deviation when applicable. 

c Values are calculated within the Jh > 36 mA m — 2 con- 
tour; the centroid is near 3 Mm in height. 

called "magnetic tongue" (e.g. iLuoni et al.l [20111 ). Fast 
sunspot rotation and flux cancellation are also observed. 

AIA 171 A observation of the extreme-ultraviolet 
(EUV) light emitting plasma roughly outlines the coro- 
nal magnetic field (Fig. QJ). As a validation of our field 
extrapolation, we plot selected field lines for comparison 
and color-code them according to \J Z \ at their footpoints 
(Fig. [TJl). Brighter color, such as red, corresponds to 
stronger \J Z \. The modeled field lines with strong cur- 
rents indicate non-potential structures, and are qualita- 
tively in good agreement with observation. In particular, 
field lines rooted in the core region (regions near the ma- 
jor PIL) morphologically resemble the observed EUV fil- 
ament. In contrast, potential-like loops (with weaker cur- 
rent) further away from the center are not well recovered, 
especially on the north side. A more detailed discussion 
on the extrapolation can be found in Appendix [Aj 

A closer look at the core region, summarized in Fig- 
ure [21 reveals a close match between the shape of the 
filament and the photospheric PIL. There seem to be a 
few faint strands in the dark filament: NLFFF extrapo- 
lation suggests that an ensemble of highly sheared loops 
thread the plasma. These loops are rooted in PI /NO in- 
side a narrow strip of strong-current, high-a distribution 
along the PIL. These loops are typically low- lying, with 
apexes well below 10 Mm. We note the broad distri- 
bution of photospheric a, mainly from to 0.5 Mm -1 
(Fig. (2]i), necessitating the non-linear treatment of any 
realistic modeling attempts. 

Figure [2^ shows in the background the computed hor- 
izontal current density (Jh) in a vertical plane perpen- 
dicular to the PIL (near "A" in panel c). The projected 
magnetic field vectors on this cross section rotate around 
an axis at about 3 Mm altitude; Jh peaks at the same 
height. Signatures of twisted flux ropes, such as an X- 
point or a hollow core Jh distribution around its axis 
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Fig. 2. Close-up view of the AR core field structure for the 2011 February 14 20:35:20 UT frame, (a) Unsharp masked AIA 171 A 
image. The FOV is denoted by the dotted box in Figure QJ;. Two straight lines show the baselines of the vertical cross sections (CS1 and 
CS2) through the computation domain plotted in panels (e) and (f). Yellow dotted lines are manually drawn to outline the dark iS'-shaped 
filament plasma. One of the faint strands appearing in the filament as well as a few bright overlying loops are marked on the image, (b) 
Selected NLFFF field lines. The torsional parameter, a, derived from preprocessed magnetogram is shown as background. Contours are 
for B z = ±600 G; the cyan line is for the PIL. (c) Preprocessed, remapped HMI vector magnetogram showing the emulated chromospheric 

field. The yellow contours indicate the 90% level of the magnetic connectivity gradient metric, N (equation IB2I in Appendix [bJ . The 
contour is an indicator of the intersection of quasi-separatrix layer (QSL) with the lower boundary. A, B, and C mark features of interest, 
(d) Histogram of a distribution in the FOV of panel (b). The peak near 0.46 Mm -1 comes from the region along the PIL. (e) Horizontal 
current density ( J^) distribution on a vertical cross section (CS1) with projected NLFFF field vectors. Only the component perpendicular 
to the cross section is shown for J^. The white and black contours are for 10 and 36 mAm -2 J^, respectively, (f) Similar to (e), showing 
inverse-polarity configuration (CS2), where field vectors turn upward from right to left. 



fe.g. lBobra et aLll2Q08t fSu et al.ll2011l ) are not obvious. 

Using AIA 304 A images (Fig. [3j) where the filament 
plasma is best defined, we estimate the width of the 
EUV filament to be about 6 Mm, slightly wider than 
it appears in the 171 A band. The apparent length is 
at least 90 Mm. We choose the 10 mAm -2 con- 
tour (Fig. Et) as a proxy for the filament cross section 
as it encloses a region comparable to the observed fil- 
ament width. A summary of the estimated field pa- 
rameters in this region (filament channel) are listed in 
Table [TJ In particular, the field strength (B) near 
the Jh maximum is about 1200 G. This inferred kilo- 



gauss field is strong com pared with previously mod- 
eled plage filaments (e.g. JAulanier fc Demoulinl [2003; 
IGuo et al. 2008; Jing et al. 2010|) and recent observation 
(e.g. Kuckein et al. 2009) . We estimate the field twist 
angle as (j) = aL/2 (jLongcope et al.lll998[ ) with L being 
the loop length and find an average of 0.6 turn, with 
about 0.9 turn near the Jh maximum. 

The field line mapping appears to "bifurcate" near the 
PIL on the east side ("B" in Fig. [2fc). Some EUV loops 
originate from the PI magnetic tongue connect back to 
NO; others deviate from the PIL and connect to Nl more 
than 50 Mm away (Fig. and [2^). The photospheric a 
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changes sign here as well. We exploit the extrapolation 
and calculate the magnetic connectivity gradient metric 
N (equation IB2I in Appendix [B]) on the lower bound- 
ary (jDemoulin et al.lll996[ ). High-iV contours show the 
photospheric intersection of the quasi-separatrix layers 
(QSLs). They divide the magnetic tongue into two dis- 
tinctive regions (Fig. where the computed field con- 
nectivities diverge. 

In additio n, an interesting; "i nverse-polarity" configu- 
ration (e.g. iMackav et al.l I2010L and references therein) 
exists nearby too, where horizontal field on the photo- 
sphere points from negative vertical polarity side to the 
positive (near "C" in Fig. Et). On the vertical cross sec- 
tion shown in Figure [2]f, projected field vectors form a 
"concave-up" configuration. Loops turn upward here or 
become nearly parallel to the photosphere. 

3.2. Evolution of Magnetic Field 

We show five representative snapshots of the evolv- 
ing AR in Figure [3] (see online animation for the en- 
tire 5-day evolution with full cadence). The left and 
middle columns show B z and negative AIA 304 A im- 
ages, respectively. The right column shows the verti- 
cal integration of absolute current density (J) over the 
lowest 10 Mm as derived from the NLFFF extrapola- 
tion. Patterns of strong current concentration may serve 
as a proxy for the n on-potential coronal structures (e.g. 
ISchriiver et alMOQl . 

Here we highlight a few key features of the magnetic 
field evolution. For reference, the total unsigned flux |<£| 
and its change rate (d|$|/dt) are plotted in Figure 2k; 
the photospheric unsigned current (|/|) in Figure 2t>; the 
GOES SXR flux in Figure 2k- 

— The early stage of the AR development (top two 
rows of Fig. [3]) features fast flux emergence that 
coincides with the appearance of a pronounced fil- 
ament. Such fast emergence continues for about 1 
day (Fig. 2k) with a rate of several 10 20 Mxh -1 , 
mainly in th e aforementioned tw o bipoles P0/N0 
and Pl/Nl (ISchriiver et al.ll20TTh . Note that NO 
consists of two sunspots, the western one emerges 
much later. AIA images show frequent brighten- 
ing as PI and NO converge and collide; new mag- 
netic connectivities are being established rapidly 
in the corona. A filament becomes visible along 
the PIL within a few hours. Along the PIL, signifi- 
cant electric current injection appears to take place 
(Figs. [3)3 andHb). During this interval, the net (un- 
balanced) current in each magnetic polarity also in- 
creases drastically, hinting t hat the newly e merged 
flux is highly non-potential (jSchrijverll2009[ ). 

— After the initial fast flux emergence (third row of 
Fig. [3]), the sunspot complex develops further with 
slower flux emergence but lasting strong shearing 
motion between PI and NO. The injection of cur- 
rent through the photosphere clearly slows down 
(Fig. 2fc>)- However, the vertically integrated \J\ 
still increases by about 25% over a period of 40 
h, suggesting a net build-up process in the corona. 

5 http : //sun . Stanford . edu/~ xudong/Article/evo . mp4 



A pair of newly emerged pores (P2/N2 in Fig. [3b) 
with a few 10 20 Mx flux undergo fast shearing on 
the northeastern side and substantially reconfigure 
the coronal field. Smaller recurrent eruptions take 
place above these pores throughout the next two 
days. At the same time, the filament appears to 
be stretched. It increases in length, becomes some- 
what warped (perhaps elevated) on the western end 
and extends further west to P0 (Fig. [3h). Simul- 
taneous flux cancellation is observed between PI 
and N3, which may act favorably to the subsequent 
eruptions. 

— Toward the end of the 5-day period (bottom two 
rows of Fig. [3]), flux emergence becomes episodic, 
sometimes overtaken by flux cancellation (Fig. 2k)- 
Nevertheless, the shearing between NO and PI con- 
tinues. The integrated \J\ in the AR core region 
shows little sign of decreasing, and is fluctuating 
within 10% of the pre-flare value. A similar trend 
is present in |/| as well (Fig. 2t>)- Over time, the 
converging motion between sunspots of the same 
polarity starts to simplify the region toward a more 
bipolar structure. The filament seems to survive 
multiple eruptions and is still visible at the end of 
the 5 days. 

3.3. Evolution of Magnetic Energy 

We now consider the distribution of the magnetic free 
energy based on the NLFFF extrapolation. Spatially, 
the low-lying, current-carrying core field demonstrates 
strong concentration of free energy in the AR core, from 
the chromosphere to the lower corona. The free energy 
density (ef) appears largely co-spatial with the current 
distribution (Fig. [3] right column). At T = h, the 
20% level contour of the peak e/ (vertically integrated) 
accounts for only 9% of the AR area (Fig.[3h), but 77% of 
the overall free energy. The height profile of €f plotted in 
FigureHf shows that about 50% of all free energy is stored 
below 6 Mm, and 75% below 11 Mm. We note that the 
maximum free energy density occurss at about 3 Mm 
altitude, corresponding to the height of the peak in Jh 
(Fig. [2k). These distribution patterns remain relatively 
stable after the filament becomes well defined in AIA 
images around -36 h. 

The temporal profiles of the total NLFFF and PF mag- 
netic energy (En and Ep) in Figure 2fc roughly scale with 
the unsigned flux. On the other hand, the estimated 
magnetic free energy Ef in Figure 2H shows interesting 
variations in time and appears sensitive to AR activity, 
as one would expect. We plot the ratio En /Ep in Fig- 
ure 2k as an additional measurement of the AR's non- 
potentiality. Uncertainties reported in this section only 
represent the effect of spectropolarimetric noise and are 
obta ined from a pseudo Monte-Carlo method (see Sec- 
tion 15. ip . Systematic uncertainties from the extrapola- 
tion algorithm may be greater. Again, several key fea- 
tures are highlighted below. 

— Rapid free- energy injection during flux emergence. 
Ef increases drastically following the flux emer- 
gence (Fig. 2k and b). Free energy concentration is 
seen in the vicinity of the filament (Fig. (3b), and e/ 



Sun et al. 



120 




02.12 19:59:22 UT 



4000 ■ 

o" 





[DN] 

AIA 304 A T = -52 h 



Filament 



02.13^07:59:22 UT 




T = -40 h 



02.14_23:59:?2 UT 




T = h 



02.15 11:59:22 UT 



T = 12h 



• 1 ■ 

T ■ 



02.1 6_1 1:59:22 UT 



T = 36 h 



02.12 19:59:19 UT 



Integrated 
Current Density 



2.6H' 

0.1 U- 

[x10 5 Am" 1 } 
T = -52h 



Filament 



02.13 07:59:20 UT 



T = -40 h 



02.14 23:59:20 UT 



T = Oh 



02.15 11:59:20 UT 




T= 12 h 



02.16 11:59:20 UT 




T = 36h 



40 80 120 160 40 80 120 160 
X (arcsec) X (arcsec) 



40 80 120 160 

X (arcsec) 



Fig. 3. — Five snapshots of the evolving AR11158. They are taken at about T = -52 h, -40 h, h, 12 h, and 36 h, with February 15 00:00 
UT as time 0. Left column: HMI B z as in native coordinate (as recorded by camera). Middle column: negative AIA 304 A image showing 
chromosphere and transition region structures in which the AR filament is best discernible. Right column: vertically integrated current 
density from NLFFF extrapolation over the lowest 10 Mm. The thick solid, thin solid, and dotted contours are for similarly integrated free 
energy density at 80%, 50%, and 20% of the peak value for frame T = h. Images are remapped back to the native coordinate for direct 
comparison of HMI and AIA observations. The box in (b) indicates the field-of-view of Figure [10] for a flux-emerging region. Features of 
interest are marked in some panels; see text for details. (An animation of the entire 5-day period is available in the online journal.) 



at all altitudes increases significantly (Fig. Hf). In 
particular, Ef shows a twelve-fold increase from -55 
h to -45 h, reaching 1.07 x 10 32 erg and amounting 
to 43% of the final maximum Ef two days later. 
The ratio of Ejy/Ep peaks at about 1.57 at -45 
h. Significant changes in the coronal structures en- 
sue, as conveyed by the AIA images (animation for 
Fig. [3]). Flux change rate and Ef then plateau for 
a few hours, and Ejy/Ep drops back to about 1.30. 



Careful inspections of the vector magnetogram re- 
veals signatures of emerging flux tubes, which we 
discuss in Section I5TT1 

Gradual energy build-up and decay. After the ini- 
tial fast increase, Ef accumulates at a slower rate 
(Fig. [Hi). The AR has attained a well-developed 
quadrupolar topology. The growth rate of Ef and 
the ratio Ejy/Ep are both nearly constant. Simi- 
lar trends can be found for the post X- flare stage, 
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Fig. 4. — Evolution of magnetic energy and related quantities of AR 11158 over 5 days, (a) Total unsigned magnetic flux (\&\) and 6-h 
smoothed flux change rate (d|<I>|/d£). (b) Total unsigned current (|/|). (c) Magnetic energy derived from the NLFFF and PF extrapolation 
(En and Ep). (d) Estimated magnetic free energy (Ef). (e) Ratio between the NLFFF and PF energy (Ejy/Ep). (f) Time-altitude 
diagram of average magnetic free energy density in the AR center (FOV of Fig. (5^f))- Dotted lines indicate the height below which the 
accumulated values reach 50% and 75% of the total. Curve on the left shows the height profile for T = h. (g) GOES 5 m soft-X ray flux 
(1-8 A channel). The 5 frames in Figure [3] are marked as circles. Flux and current are derived using only pixels with \B Z \ greater than 
100 G. Error bars in (c)-(e) show the estimated effect of noise and are plotted every 6 h, but are usually too small to be seen. Insets of (d) 
and (e) show results with a 12 minute cadence from -1 h to 8 h (shaded grey band). Error bars right before and after the X-class flare in 
these insets are more visible. The vertical dotted lines indicate the peak time of the X-class flare at 1.9 h. 
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except now Ef starts to decrease and the region 
relaxes to a more potential state (Fig. 0]i and e). 
The energy dissipation rate in the corona probably 
has exceeded the growth rate from the photosphere 
as the flux emergence slows down (Fig. |4^l). 

— Step-wise energy loss during the X-class flare. Ef 
reaches (2.47 ±0.03) x 10 32 erg prior to the X-class 
flare, accounting for about 23% of the total energy. 
During the flare, Ef displays a step-wise sudden 
decrease (Fig. [Hi and inset). The 1 h average of Ef 
before and after the flare differ by about (0.34 db 
0.04) x 10 32 erg, 14% of the pre-flare Ef. This 
value situates at the lower end of what is adequate 
to power a typical X-class flare (e.g. [Hudson 20111 
and references therein) and may be intrinsically an 
underestimation, as discussed in Section 15.11 A 
clear discontinuity also appears in the Ejy/Ep and 
€f height profile. Sudden energy decrease is also 
found during the earlier M-6.6 flare (-30.4 h) with 
a smaller magnitude. 

The continuous monitoring of the AR free energy with 
relatively high temporal/spatial resolution is made pos- 
sible, for the first time, by the HMI vector field obser- 
vations. This is especially useful for the study of ma- 
jor eruptive events. Due to the nature of the extrap- 
olation method, changes in Ef here are determined by 
the boundary conditions. Therefore, the step-wise en- 
ergy loss found during the X-class flare must be related 
to rapid and significant field changes on the photosphere. 
We consider this change in detail in Section [U 

4. FLARE-RELATED RAPID CHANGE 

Figure [5] illustrates the rapid magnetic field changes 
during the X-class flare. Strong and perman ent enhance- 
ment of Bh in the AR core was reported by I Wang et~a71 
(|2011[ ) based on HMI linear polarization signal, and is 
confirmed here by full Stokes inversion. From 01:35 UT 
to 02:11 UT, the mean B h along PIL (Fig. [5fc and d) 
increases from about 1200 G to over 1500 G, 28% within 
0.6 h. Changes in B z appear less significant (Fig.[5t and 
f): mean \B Z \ decreases by about 5%. We skip the two 
frames in between (01:47 and 01:59 UT) t o av oid possible 
artifacts from flare emission (see Section [572]) . 

Shown in Figure El the difference image of pre- 
and post-flare Bh bears a striking resemblance to the 
Ha flare rib bons observed by th e Solar Optical Tele- 
scop e (SOT; Tsuneta et all [2008) on the Hinode satel- 
lite (|Kosugi et al.ll2QQ7h . The elongated regions that are 
swept by the evolving ribbons or lying in between show 
significant Bh enhancement, whereas patches with de- 
cayed Bh appear on both the north and south sides. In 
the narrow region along the PIL, distribution of Bh shifts 
by about 300 G in the histogram (Fig. [7h), displaying a 
strong boost in the kilogauss range. On the other hand, 
\B Z \ often decreases where Bh increases, but the signals 
are weaker and appear to be mixed with the opposite. 
The distribution of \B Z \ gently decays in the strongest 
part, but otherwise remains similar (Fig. 01). 

The combined effect is that the field becomes overall 
stronger and more inclined in the AR core (Fig. O). In 



TABLE 2 

Flare-related change in magnetic energy and photospheric 

FIELD. 



Parameter* 1 


Unit 


Pre-flare b 


Post-flare b 


Difference 


Ef 


10 32 erg 


2.47 ±0.03 


2.13 ±0.03 


-0.34 ±0.04 


En/Ep 




1.29 ±0.00 


1.25 ±0.00 


-0.04 ±0.01 


(B h ) 


10 3 G 


1.20 ±0.02 


1.53 ±0.01 


±28% 


(\Bz\) 


10 3 G 


0.96 ±0.03 


0.91 ±0.03 


-5% 


(B) 


10 3 G 


1.66 ±0.01 


1.87 ±0.01 


±13% 


(M) 


degree 


36 ±1 


29 ± 1 


-19% 


(\J Z \) 


mA m -2 


36 ±2 


27 ± 1 


-25% 


(a) 


Mm- 1 


0.48 ± 0.03 


0.30 ±0.02 


-37% 



a Ef and the ratio En/Ep are computed using the center region 
(184 x 144 x 115 Mm 3 ) of the extrapolation domain. Other param- 
eters are for the AR core photospheric field within the boxed region 
in Figure [5J:. Notations are same as Table [T] 

b For Ef and E^/Ep, pre-flare values are calculated from the av- 
erage of 5 frames between 00:47 to 01:35; post-flare between 02:11 
and 02:59. Uncertainties are estimated effect from noise. For all oth- 
ers, pre-flare values are calculated for the 01:35 frame; post-flare for 
02:11. Values are reported as mean ± standard error for comparison. 
c Difference for Ef or Ejy/Ep is the absolute value, others are per- 
centage differences. 

addition, the azimuth of Bh appear to change in a fashion 
such that they become better aligned and more parallel 
to the PIL in its vicin ity (Fig. [T^t a nd b), consistent with 
previous reports (e.g. Wa ng et al.lll994l ). After the flare, 
a appears to be smaller, suggesting the field is perhaps 
less twisted. A summary of the pre- and post-flare field 
parameters is provided in Table [2 

These fast changes suggest a scenario in which the 
change of coronal connectivity, probably induced by re- 
connection, feeds back to the photosphere. In particular, 
newly reconnected loops with footpoints located in the 
flare ribbons, both shorter and more parallel to the PIL, 
are a priori consistent with t he ob servations. We discuss 
the topic further in Section 15.21 Morphologically, the 
changes are in line wit h the conjectured magnetic "im- 
plosion" (Hudson 2000): a decrease in coronal magnetic 
energy during explosive events should lead the coronal 
field to contrac t, resulting in a "m ore horizontal" pho- 
tospheric field (|Hudson et al. I [20081 ). Here we analyze a 
series of NLFFF snapshots and monitor the Jh distribu- 
tion on a vertical cross section through the computation 
domain (Fig. [5^ and b). The patterns show an apparent 
contracting motion at flaring time owing to the altered 
boundary condition (see online animation). The apex 
of the outermost Jh contour lowers by about 3 Mm in 
Figure Eb. 

Motions of coronal loops farther away from the AR core 
provide evidence for the conjecture. During flares and 
CMEs, a depletion of magnetic energy leads to smaller 
magnetic pressure gradie nt, and loops must contract to 
reach a new balance (e.g. [Liu fc Wangll2010[ ). We show 
AIA 171 A observation for this event in Figure El Ex- 
pansion of the coronal structure and faint circular propa- 
gating fronts are visible startin g around 01:48 UT, p rob- 
ably linked to CME initiation (jSchriiver et al.ll2Qll ). At 
about 01:50 UT, the southern potential-like loops sud- 
denly move inward; the northern loops follow immedi- 
ately (see online animation for Figure [8]). The retraction 
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Fig. 5. — Rapid field changes at the photosphere and in the corona during the 2011 February 15 X-class flare, (a) Jh distribution on 
a vertical cross section as derived from NLFFF extrapolation, before the flare at 01:35:20 UT. Only the component perpendicular to the 
cross section is included. The location of the cross section is indicated in panels (c)-(f) as a long straight line. The dotted, dashed and solid 
contours indicate values of 2, 10, and 36 mAm -2 respectively, (b) Same as (a), for 02:11:20 UT after the flare, (c) Remapped observed 
by HMI for 01:35:20 UT as viewed from overhead. The dotted (solid) lines show contour for 1600 G (1200 G). Places with significant field 
change are marked by arrows and boxes, (d) Same as (c), for 02:11:20 UT. (e) Remapped B z at 01:35:20 UT. The dotted (solid) contours 
are for ±2000 G (±1000 G). (f) Same as (e), for 02:11:20 UT. The large yellow box indicates the region evaluated in Figure EJd) . (An 
animation of this figure is available in the online journal.) 
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Fig. 6. — Difference images of pre- and post-flare magnetic field, (a) Difference of Bh between 02:11:20 UT and 01:35:20 UT. The two 
frames are kept in the native coordinates without remapping, but are co-aligned with Carrington rotation rate to the nearest full pixel. 
Contours are for the mean B z of the two frames, at ±600 G. (b) Difference image of \B Z \. (c) Ha flare ribbons observed by Hinode /SOT 
at 01:53:14 UT, near flare peak and roughly midway between the two HMI frames. 
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Fig. 7. — Close-up view of the pre- and post-flare vector field in the AR core region, (a) Remapped vector magnetogram for 01:35:20 
UT. The yellow contour is the PIL; the boxed region is identical to Fig. 0;. (b) Same as (a), for 02:11:20 UT. (c) Histograms of at two 
different times, in the boxed region in panel (a), (d) Histograms for \B Z \. (e) Histograms for inclination angle | . 
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Fig. 8. — Observed EUV coronal loop retractions, (a) AIA 171 A image at 01:48:01 UT on February 15, after the onset of flare. Two slits 
that are largely perpendicular to the loops are used to obtain the time-position diagrams in the following panels. Arrows along the slits 
indicate the approximate direction of the transverse motion, (b) Time-position diagram for slit 1 constructed by stacking a time sequence 
of coaligned image slices from left to right. Only the images with the normal exposure time are used, resulting in a 24 s cadence. Images 
are coaligned to sub-pixel accuracy. The dotted line shows an inward loop contraction pattern with a transverse speed of 43 kms -1 . The 
horizontal patterns in the upper half are from the features in the background. Loop oscillations are visible, (c) Same as (b), for slit 2. 
Dotted line shows a transverse speed of 18 kms -1 . Patterns of expanding loops also appear around 01:50 UT, moving from position 20 
Mm toward 0. (An animation of this figure is available in the online journal.) 
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Fig. 9. — Altitude profile for extrapolated field around the time of the X-class flare, (a) Profiles of average magnetic free energy density 
in the AR core region (FOV of Figure [5f) at February 15 01:35:20 UT and 02:11:20 UT. Shaded grey band indicates the uncertainty due 
noise, and is usually to narrow to be seen, (b) Time-altitude diagram for free energy density during the 9 h around the X-class flare, similar 
to Fig. [4f. The vertical dashed line indicates the flare peak time at 1.9 h. Green and pink vertical lines indicate the two profiles shown in 
panel (a), (c) Profiles of weighted mean shear (6 W ) near the PIL (boxed region in Figure [5f), at February 15 01:35:20 UT and 02:11:20 UT. 
(d) Time-altitude diagram for W . Dotted lines show the heights at which the shear reaches 15°, 30°, and 45° respectively. PF becomes 
mostly east-west directed at about 20 Mm, a height much lower than that of NLFFF, thus resulting in the 6 W increase above 20 Mm. 



proceeds for about 5 minutes, mostly during the flare im- 
pulsive phase. We place two slits near the loop apexes to 
construct time-position diagrams (Fig. [8b and c). The 
inward motion shows a projected speed on the order of 
tens of kilometers per second, some with apparent oscil- 
lations. The transverse displacement can be as large as 
15 Mm. The possible ambiguities in these observations 
are discussed and resolved in Section [5^21 

According to NLFFF extrapolation, such contraction 
yields a more "compact" energy distribution and a coro- 
nal field that relaxes more rapidly with height. We show 
in Figure [9^ two altitude profiles of the mean free en- 
ergy density before and after the flare. While the total 
magnetic free energy decreases after the flare, a larger 
percentage gets stored in the lower altitudes. Higher 
overlying field thus becomes even less energetic. Simi- 
lar trends can be found in the field- weighted mean shear 
W (equation IB4I in Appendix [BJ. This quantity mea- 
sures the field- weighted direction difference between the 
NLFFF and PF field, thus providing a quantitative de- 
scription of the non-potentiality of the magnetic field. As 
shown in Figure [9k, the whole profile moves downwards 
by about 2 Mm, equivalent to a shear decrease at almost 
all heights. The bottom layer, as an exception, becomes 
slightly more s heared because B h is now more parallel 
to the PIL (cf. iWang et!Hll994f ). The time-height dia- 
grams covering 9 h around the flare with 12 m cadence 
more clearly convey the sudden and permanent nature of 
the changes (Fig. and d). 

5. DISCUSSION 



5.1. On the Coronal Energy Budget 

As discussed in Section O early flux emergence brings 
along a large amount of electric current and free energy. 
We show one of the flux emergence sites in Figure [TO1 
Highly inclined fields (inclination angle 7 measured from 
the photosphere less than 20°) appear within a widen- 
ing patch, separating two vertical flux concentrations of 
opposite polarity that themselves grow in size and field 
strength in time. These patterns appears consistent with 
an emerging flux-tube at its initial stage of rising into the 
solar atmosphere. 

Before -45 h (February 13 03:00 UT), E N increases 
rapidly with the current injection from the increasing 
Bh. These highly inclined fields do not cause vertical 
flux to increase proportionally. Since the extrapolated 
potential field solely depends on B Z1 the growth of Ep 
lags behind. This may explain the early peak in En /Ep 
shown in Figure 2k . The decrease of En /Ep can in turn 
be explained by the subsequent fast increase in vertical 
flux relative to the horizontal counterpart when the axis 
of the flux tube approaches the photosphere. 

We estimate the effect of spectropolarimetric noise on 
the energy estimation by conducting a pseudo Monte- 
Carlo experiment. In general, the noise level for the 
inverted HMI longitudinal field is below 10 G; it is 
about an order highe r for the transverse component 
(Hoekse ma et al.ll2012[ ). For each frame, we create 10 
extra-noisy Stokes profiles by adding to them a HMI- 
characteristic level of photon noise (on the order of 0.1% 
of quiet Sun intensity). Inversion, disambiguation and 
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Fig. 10. — Three snapshots of a flux emergence site. The snapshots are 3 h apart. The FOV corresponds to the boxed region in FigureEJ). 
B z is plotted in the background; contours are for ±600 G. Arrows indicate highly inclined B^. For clarity, arrows are drawn only on every 
third pixel, for field strength 200 G < B < 800 G, and for inclination \^\ < 20°. Images are in native coordinates. 



modeling are subsequently performed on these 10 inputs, 
and the standard deviation is adopted. The estimated 
noise-related uncertainty in En is typically around 0.6%, 
accounting for 2%-4% of Ef (Fig. |4]). This result is 
consistent with earlier tests on analytical field configu- 
rations and synth etic spectra with larger run numbers 
(T iwari et a"Lll2009[ ). The relatively small effect is proba- 
bly owing to the preprocessing procedure, which ensures 
force-freeness by smoothing away the spurious compo- 
nents. We nevertheless caution that the percentage un- 
certainty may be above 10% for Ef during the early stage 
of AR development, when the field is weaker and the 
measurement signal-to- noise (S/N) is lower. 

To corroborate the result from the NLFFF extrapola- 
tion, we apply the magnetic vi rial theorem to the pre- 
processed lower boundary (e.g. Met calf et al. 2005, and 
references therein; see also equation IB5I in Appendix IB]) . 
If the boundary is force-free enough and the region is 
largely isolated, as is the case for this HMI time series 
(see Appendix |A|) , the free energy in the volu me is ex- 
pected to be well recovered ([Metcalf et al. 2008). Around 
the X-class flare, the evolution of En computed from the 
virial theorem shows a similar trend to that from the 
NLFFF extrapolation, as shown in Figure [TTk : the en- 
ergy decreases during the flare are similar as well. En 
derived from the virial theorem is about 14% higher. 

For the energetics during the flare, we independently 
estimate the nonthermal electron energy using the hard 
X-ray (HXR) spectra from the Ram aty High Energy So- 
lar Spectroscopic Imager (RHESSI; iLin et al. 2002 ) un- 
der the classical thick-target model (e.g. iHolmann et al.l 
I2003L and references therein). The nonthermal electron 
flux spectra can be well fitted by a very steep power law 
with the index 5^5-10 and a low energy cutoff £? c ~20-30 
keV. The total nonthermal energy integrated from 01:44 
UT to 02:17 UT is estimated to be -0.54 x 10 32 erg, 
as shown in Figure HTb . Since the energy needed for 
thermal plasma heating, CME lift-off, and flare radia- 
tion are expected to be on the same order of magni- 
tude, if not greater than the nonthermal electron energy 
(jEmslie et al.ll2004t ). the estimated total energy budget 
has to be on the order of 10 32 erg. 

Based to this analysis, the computed pre-flare Ef 
(2.47 x 10 32 erg) is actually adequate to power this 
flare and accompanying CME eruption. On the other 
hand, the difference between the pre- and post-flare Ef 
(0.34 x 10 32 erg), which we used as the estimate for en- 
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Fig. 11. — Estimation of flare energy release, (a) Magnetic en- 
ergy estimated from the NLFFF extrapolation (same as Fig. [4j 
and the magnetic virial theorem. The magnetic virial theorem is 
applied on preprocessed lower boundary; the uncertainty is derived 
from a Monte-Carlo approach. Most error bars are too small to be 
seen. The grey background indicates the period during which the 
nonthermal electron energy is fitted in the next panel, (b) Accu- 
mulated nonthermal electron energy derived from RHESSI HXR 
spectra using the thick-target model. Spectra for seven individual 
time intervals covering 01:44 UT to 02:17 UT (grey background) 
are fitted, the time range of each is noted by the horizontal bars. 
The HRX light curve in the 25-50 keV band is plotted in the back- 
ground on a logarithmic scale. The vertical dotted line (01:56 UT) 
shows the SXR peak in GOES light curve. 



ergy loss, appears too small. One speculation is that 
the post-eruption coronal field is extremely dynamic and 
possibly deviates from the force-free condition, thus can- 
not be well described by a static force- free snapshot. If 
the field is locally (immediately above AR core) less en- 
ergetic than what the NLFFF extrapolation indicates, 
then the discrepancy in the energy estimation will be al- 
leviated. Another speculation involves the preprocessing 
scheme where a certain amount of smoothing is applied 
to the photospheric data to ensure the force-freeness of 
the boundary. Smoothing generally reduces the elec- 
tric current inferred from the boundary, thus tending to 
lower the energy in the volume. Coordinated chromo- 
spheric magnetic field measurement from other obser- 
vatories may help to constrain this procedure. A third 
possibility involves the photospheric field measurement 
itself. Small-scale field structures (sub-arc second) with 
high field gradients may carry strong electric current and 
a non-negligible amount of free energy, but will be unre- 
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solved with the moderate spatial resolution of HMI (1"), 
and thus overlooked in the modeling result. 

We note that a series of comparative and diagnostic 
efforts pointed out the difficulties in em ploying realistic 
data input for NLFFF ex trapolation (e.g lSchrijver et all 
2008; DeRosa et al. 2009). Different algorithms can pro- 
duce very different energy estimations and field topolo- 
gies, which mainly arises from the inconsistency between 
the input boundary condition and the force-free assump- 
tion (Appendix [Aj). We point out that the selection of 
the boundary data, e.g. the size of FOV, the degree of 
flux balance, can a lso directly affect th e quality of the ex- 
trapolation re sult (DeRosa et al. 2009) . Recent works by 
IWiegelmann fc Inhesterl (|2Q10h and IWheatland fc Lekal 
(|2011f ) included observational uncertainty in the revised 
algorithms. The forme r was tested with HMI data 
(Wiegelmann et al. 2011) and the later with SOLIS data 
(Gilchrist et al. 2011). Both have shown to improve the 
quality /consistency of the solution. Further advance- 
ment of the algorithm is necessary and seems promising. 
Nevertheless, due to the aforementioned discrepancies, 
we think that quantitative energy estimations derived 
from modeling should be interpreted with care in the 
meanwhile. 

5.2. On the Flare- Related Field Changes 

The rapid (~10 min) field change at the flare time 
likely has a coronal driver, as the photo spheric Alfve n 
timescale is too long to be compatible (e.g. lHudsonll201 lb . 
Reconnection, in particular, may change the coronal field 
connectivity fast enough to result in large-scale field vari- 
ation at the "line-tied" loop footpoints. We illustrate one 
of the consequences in Figure [12] using the NLFFF ex- 
trapolation. Before the flare, highly sheared field lines 
(Fig. fT2k ) cross over each other near their footpoints, 
creating a configuration in favor of tether-cutting recon- 
nection. Sigmoidal core loops "are sheared past each 
other so that the y overlap and are crossed low above 
the neutral line" (jMoore et al.ll200T| ). These loops ex- 
change footpoints during the reconnection, and the inner 
footpoints indeed appear to be well linked in the post- 
flare field in the NLFFF extrapolation (Fig. fT2b). New 
loops become shorter, consistent with the enhanced Bh 
on the boundary. We note that the illustration is purely 
schematic without any reference to the reconnection pro- 
cess itself. Moreover, loops on the north and south sides, 
further away from the PIL, may become connected after 
the flare, resulting in a more vertical configuration where 
Bh decreases (Fig. EH)? simi lar to the e f fect o f £-spot um- 
bral darkening reported in iLiu et al.l (j2QQ5h . This sce- 
nario has been indirectly probed thr ough the contrac- 
tion of HXR sources (jJi fc Wa ng 2007) a nd the mishear- 
ing m otion of EUV flare ribbons (|Ji et al.ll2006l : ISu et al.l 
120071 ), Direct imaging can be difficult as EUV or SXR 
images are usually saturated during major flares. The 
complicated temperature structures of the flaring atmo- 
sphere may also prevent a straightforward interpretation 
of the observations. 

Besides the short loops near the PIL, tether-cutting 
reconnection also produces a flux-rope that connects the 
two far ends of the sigmoid after the exchange of the foot- 
points (jMoore et al. 2001). In a detailed account of the 
coronal activitie s in this event, Schriive r et al.l (j2011f ) use 
MHD modeling (jAulanier et al .11201 Of ) to interpret the ob- 



servations and find that the flux-rope formed by recon- 
nection eventually erupts as part of the CME structure. 
Interesting evidence a lso arises from the "s unquake" ob- 
servation report ed by Zharkov et al.l (|2Qllh following the 
initial report by lKosov ichev (2 01 If ). Two distinctive pho- 
tospheric seismic sources are found at the two far ends 
of the sigmoid and their initiations closely match the 
flare/CME onset time. Because the signals appear before 
the HXR emission peak time and do not spatially match 
the HXR sources, they cannot be explained by particle 
precipitation as many earlier events are. Instead, flux- 
rope eruption appears to be the trigger, and magnetic 
variation is proposed as one of the possible mechanisms. 
The coronal fields rapidly reconstruct after the eruption 
and the overall effect may ge nerate a strong dow nward 
impulse into the photosphere (Hu dson et al.ll2008h . Nev- 
ertheless, it is pointed out that the seismic sources do not 
appear to be co-spatial with the strongest magneti c field 
or the strongest field variations (|Kosovichevll2011f ). The 
exact excitation mechanism remains to be explained. 

Regarding the possible ambiguity in the loop retrac- 
tion observed in AIA 171 A (Fig. [8j), we can rule out 
the scenario where loops remain steady but only cool 
sequentially into the EUV temperature sensitive pass- 
band. Such an explanation cannot account for the sud- 
den start and stop of the motion, or the apparent loop 
oscillations. To rule out the projection effect, we inspect 
the coronal image s equences taken by SECCHI EUVI 
([Howard et al.l 12008( 1 on the twin STEREO spacecraft 
that are both near quadrature from the Sun-Earth line. 
Coronal loops, when viewed from the side, appears to 
become "flattened" and in generally pressed down to- 
ward the solar surface during the eruption. The appar- 
ent loop length, if unchanged, should increase instead. 
Therefore the retraction is likely to be real, its ampli- 
tude and speed underestimated from overhead AIA ob- 
servation. Detailed analysis of EUV loop motion may 
provide quantitative diagnosti cs of the coronal field (e.g 
lAschwanden fc Schrijverll2011f h but is out of the scope of 
this work. 

We note that the field change discussed here must be 
distinguished from the transien t variations in the longi- 
tudinal field (Ko sovichevl l201lf ). The latter lasts only 
a couple of tens of minutes and is likely an artifact 
caused by flare-related emissions. HMI spectra at these 
abnormal pixels briefly deviate from the usual absorp- 
tion profile, showing an enhanced line center as expected 
(jMaurva et al.l 1201 If ). 

6. SUMMARY 

We have studied the magnetic field and its energy of 
AR 11158 over a period of 5 days using a series of HMI 
vector magnetograms. A NLFFF extrapolation, coupled 
with coronal imaging, provides information of the coro- 
nal magnetic structure, electric currents and free energy. 
From its early flux emergence to recurrent major flares 
and CME eruptions, the AR displays distinctive stages 
of energy build-up and release, driven by or resulting in 
gradual or sudden magnetic field changes. We summarize 
the major results as follows: 

— The quadrupolar AR primarily consists of two in- 
teracting bipoles. The trailing polarity of the lead- 
ing pair and the leading polarity of the trailing 
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Fig. 12. — Schematic illustration of the field connectivity as possible results from reconnection in the AR core region, (a) Selected NLFFF 
field lines prior to the flare, plotted on preprocessed vertical field. The contours are drawn at ±1000 and ±2000 G. The eastern set of loops 
(yellow) are expected to reconnect with the central set (cyan) and exchange footpoints. Side view of the loops from the left (east) is plotted 
as inset, (b) Selected post-flare loops (green) with footpoints that belonged to the two different sets of pre-flare loops in the previous panel. 
The loops become much shorter. These snapshots demonstrate that the reconnection scenario is consistent with the observed photospheric 
field. They are not intended for identifying individual reconnecting loop pairs, nor modeling the reconnection process. 



pair undergo significant shearing. A pronounced 
filament appears over the major PIL early and its 
main part persists through multiple eruptions, in- 
cluding the X-2.2 class flare and halo CME on 2011 
February 15. NLFFF extrapolation suggests that 
the field in the filament channel carries strong cur- 
rent, and is highly sheared with about 0.9 turn near 
the axis. 

NLFFF extrapolation indicates significant electric 
current and free energy injection during early flux 
emergence, about 10 32 erg over a mere 10 h. Cur- 
rent and energy mostly concentrate in or near the 
filament channel in the low atmosphere: over 50% 
is stored below 6 Mm. The computed peak free en- 
ergy reaches about 2.59 x 10 32 erg. A 0.34 x 10 32 erg 
decrease is found within 1 hour after the X-class 
flare. We show that the effect of random noise 
is small, but the systematic uncertainties can be 
large. The energy loss is most likely underesti- 
mated. 

Rapid and irreversible enhancement of the horizon- 
tal field takes place along the PIL during the X- 
class flare. The increase in is 28% on average; 
only a 5% decrease is found for \B Z \ and the signal 
is mixed. The observed photospheric field becomes 
overall stronger, more inclined and better aligned 
with the PIL. The short time scale indicates a coro- 
nal driver. Shorter loops created by tether-cutting 
type reconnection are a priori consistent with the 
photospheric signatures. The reconnection picture 
is supported by the fact that the change largely 
happens in the area that are either swept by the 
flare ribbons or lying in between. 

NLFFF extrapolation demonstrates that the 



change in the photospheric and coronal field is mor- 
phologically consistent with the "magnetic implo- 
sion" conjecture. Energy loss during the explosions 
cause the coronal field to contract to reach a new 
balance. This scenario is supported by the EUV 
loop retractions observed from AIA. The extrapo- 
lated field appears to be more "compact" after the 
flare, the lowest layer is more sheared but it relaxes 
faster with height and is overall less energetic. 

Regular quantitative study and statistical surveys of 
a large ensemble of AR vector magnetic fields now be- 
come possible with HMI observing the full solar disk at 
high spatial and temporal resolution. Although prob- 
lems still abound in both data reduction and modeling 
procedures, we anticipate more in-depth studies using 
the newly available data will eventually lead to better 
understanding of the AR magnetic fields and energetics. 
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APPENDIX 

A. EVALUATION AND DISCUSSION OF THE EXTRAPOLATION ALGORITHM 

Previous studies of coronal field extra polati on have shown the importance of appropriate input boundary conditions 
(e.g. iMetcalf et al.l 120081 : iDeRosa et al.ll2QQ9f ). It is suggested that the input magnetograms should: i) have a large 
FOV with sufficient information on the horizontal component even in the weak field region; ii) have balanced magnetic 
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flux; and iii) have balanced Lorentz force and torque. The HMI vector magnetograms used in this study appear to 
be a good candidate. For the February 14 20:35 UT frame, the region of interest with strong field (\B Z \ > 100 G) 
covers the center 11% area of the FOV (a total of 216 x 216 Mm 2 ) and is isolated from the side boundaries. There 
are no extended plage regions or ARs nearby, and the large outskirts included correspond to a coronal volume that 
contains most of the observed EUV loops. The net magnetic flux accounts for only 0.3% of the total unsigned flux for 
the strong field region, thus the flux is well balanced. In particular, the surface integrated Lorentz force and torque 
are 5-15 times lower tha n the several data sets used in earlier studies (from different instruments, for different ARs) 
(Wiegel mann et aT1l2011f ). We note that such low Lorentz force / torque may be a special case, a s the photospheric 
field is generally not expected to be force- free (Metcal f et al.|[l9 95). For the preprocessing scheme (Wiegelma nn et all 
2006) applied before the extrapolation, we test and adopt the following parameter set: \i\ = \±2 = 1, /i3 = 0.001, and 
/i4 = 0.01. Here \i\ and H2 control the net force and torque on the boundary respectively; [1% controls the relative 
influence of the observed data, and the level of smoothing. 

We evaluate the quality of the extrapolation by computing three metrics: the mean Lorentz force Lj, the mean field 
divergence L^, and the current weighted mean angle between the magnetic field and electric current ay. 

iVxtVxS), \ (A1) 



B 2 

L d = (\V-B\ 2 ), (A2) 

<T , =sin " 1 (( R ^) /(J) )' (A3) 

where B is the NLFFF field vector, B = |£?|, J = | J|, and angle brackets denote the mean value within the domain. 
Following literature, we normalize L/, to a unit volume and normalize further by 1 G 2 . Ideally all should vanish. 
For a total of 120 hourly frames computed in this study, we find Lf = 6.4 ± 2.2, = 3.1 ± 1.2, and aj = 18.0° ±1.2°. 
These values are typical for t his optimization-base d algorithm and are comparable to previously reported results based 
on other solar data sets (e.g. iSchrijver et a71l2008[ ). 

It has been pointed out that even after preproc essing, the observed field, with its uncertainties, is still incompatible 
with a strictly force- free field (jDeRosa et al.l 120091 ). The fact that a better suited input does not lead to noticeable 
improvement in this case may suggest an internal limit to the current algorithms when dealing with realistic, imperfect 
boun dary conditions. In light of this, a new implementation of the optimization scheme (Wiegel mann fc Inhesterl 
2010) has been developed and tested for HMI data. After including information on the horizontal field measure ment 
uncertainties, the solution proves to satisfy the force- free condition significantly better (j Wiegelmann et afll2011f ). For 
the 20:35 UT frame, aj is reduced from 16.6° to 5.7°. 

Problems may also arise when we adopt a large computation domain but still assume a planar boundary. A 
magnetogram spanning about 18° in Heliographic longitude, as used in this study, will then have its edge elevated by 
(sec 9° — 1)Rq (8.7 Mm, or 12 pixels) above the solar surface. This may be a reason that we fail to faithfully reproduce 
the long loops on the north side. These loops extend far away from the AR center and appear to be low- lying judged 
from the STEREO /EUVI observations. They are thus more subject to inconsistencies in the boundary condition. 
We note that the free energy estimation may be less effected in this case, as these loops are largely potential-like. 
Nevertheless, as multiple ARs are often m agnetically connected in a global scale, implementation of the model in 
spherical geometry (e.g. Tades se et al.ll2011[ ) becomes necessary. 

B. COMPUTATION OF VARIOUS MAGNETIC FIELD PROPERTIES 

The torsion parameter, a, indicates the twist and non-potentiality of the force- free magnetic structure. It may 
be computed from any component of Equation [TJ For its distribution on the photosphere, we use the vertical (z) 
component. For the filament cross section shown in Figured, assuming the normal vector of the vertical plane is n 
(which itself is horizontal), we use the horizontal component: 

^(Vx^n 

B h n y J 

The N metric ([Demoulin et al.lll996f ) is used to quantify the spatial gradient in magnetic connectivity. When the 
mapping of magnetic field line diverges, N increases drastically and indicates the existence of a quasi-separatrix layer 
(QSL). For a footpoint P(x, y, 0) on the bottom boundary in the domain, we integrate the field line passing through 
it and obtain the other footpoints (x' , y', 0). Field lines passing through the side walls or upper boundary are treated 
as missing values in this work. With the displacement vector defined as {Xi,X2} = {x f — x, y' — ?/}, N is defined as: 



N(x,y,0) 



\ i=l,2 



(B2) 



\ dx J \ dy 

The shear angle measures the angular deviation between a field vector from the potential counterpart. Let Bn and 
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Bp denote the NLFFF and PF vector, \Bn\ = I? at, \Bp\ = Bp, the mean shear angle at each point is defined as: 

-i B N • Bp 

= cos 5 5 • ( B3 ) 
The weighted mean shear # w is computed by summing over all pixels in the domain (Wa ng et al .11 19941 ): 



N 



(B4) 



The magnetic virial theorem can be used to calculate the force-free field energy in a volume solely from its boundary. 
When applied to the lower boundary {z = zo) in a Cartesian coordinate, it can be expressed as 

E = -jj- I (xB x + yB y ) B z dx dy, (B5) 

47r Jz=zo 

where E is the force-free field energy, and x, y are the coordinate of where the field is measured. Following Met calf et al.l 
(2005), we estimate the uncertainty in E through a pseudo- Monte Carlo method. We randomly vary the origin point 
of the coordinate system and repeat the evaluation for large numbers. The standard deviation is adopted. It should 
be if the field is perfectly force- free, and no magnetic flux connects outside the domain. 
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